Global existence and causality for a transmission 
problem with a repulsive nonlinearity 
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^ [ Abstract 

It is well-known that the solution of the classical linear wave equation with compactly 
supported initial condition and vanishing initial velocity is also compactly supported 
in a set depending on time : the support of the solution at time t is causally related 
<^ to that of the initially given condition. Reed and Simon have shown that for a 

real-valued Klein-Gordon equation with (nonlinear) right-hand side —Am 3 , causality 
still holds. We show the same property for a one-dimensional Klein-Gordon problem 
but with transmission and with a more general repulsive nonlinear right-hand side 
F. We also prove the global existence of a solution using the repulsiveness of F. In 
the particular case F{u) = —Am 3 , the problem is a physical model for a quantum 
>• | particle submitted to self-interaction and to a potential step. 
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1 Introduction 

For classical linear wave equations, causality is satisfied : if the initial condition has 
compact support £ and if the initial velocity vanishes, then 

supp(w(t, •)) c C(E, t) = {x e 2R/dist(x, E) < ct} 

where c is the wave velocity. The support stays in the light cone (cf. [13], Section 
X, comment on the wave equation). 

This well-known result has been generalized by Reed and Simon (cf. [28]) to the 
case of dispersive waves in three dimensions and with the nonlinear right-hand side 
—Am 3 . Our aim here is to prove the same properties for a transmission problem. The 
waves are still dispersive, the right-hand side is more general (a nonlinear function 
F) but we restrict ourselves to dimension one since we want to use the exact expres- 
sion of the solution of the linearized problem (with F = 0) which we have studied 
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in previous papers (cf. [6], [7], [9], [10]). A two-dimensional transmission problem 
is studied in [29] which could be analyzed with the methods developed below. 

Let us first introduce the mathematical context. We call Q ± — Q 2 — (0;+oo), 
I = IR + and denote by c and a k (k = 1; 2) three positive constants with a 2 > a\. 
We consider the nonlinear problem (NLP) of finding u k : I x Q k ^ IR with the time 
variable t E I and the space variable x E Q k , k — 1; 2 satisfying : 



(E k ) : ^fax) - c 2 ^(t,x) + a k u k (t,x) = F(u k (t,x)),V(t,x) E IR + * x il k 
(T ): Ul (t,0)=u 2 (t,of,VteI 

(Tx): ^(t > 0*) + ^(t > 0+) = 0,V* e J 
(/d): wjb(0, -) = /fc,A; = l;2 

(IC 2 ): ^(0,.) = 0,* = 1;2 
where F is defined from ZR to Zff and from to IR for fc G {1; 2}. 

The functions represent the initial condition and F contains the nonlinearity of 
the problem. The solution formula of the linearized problem (obtained with F = 0) 
is first obtained by Weidmann (cf. [36]) and given in [7] and [9]. The meaning of 
the coefficients is also explained in the latter two papers : the propagation velocity 
of the waves (phase velocity) is c (in the case a k = 0), and the a k s are coefficients 
which characterize the dispersion. (T ) and (Ti) express the absence of energy loss 
at the central node x = (cf. Theorem 3, Section 3 of [9]). See Section 5 of the same 
paper for a physical justification and some more details about the model. Let us just 
recall that the linearized problem is a model for a quantum particle submitted to a 
potential step. In the nonlinear problem the particular nonlinearity F(u) = —Xu 3 
means that the particle is also submitted to self-interaction as well as to a potential 
step. Cf. Hladik and Chrysos's book about quantum mechanics for an introduction 
to the physical applications of our evolution problem ([20]). 

The first part of our work is devoted to the existence and uniqueness of global 
solutions for a transmission problem with two semi-infinite axes connected at one 
point. The problem is first expressed in an abstract way (Section 2.1) and the local 
existence and smoothness of solutions is stated using the theorems stated in Reed 
and Simon's book (cf. [28], section X.13) which are based on the Banach fixed point 
theorem (cf. Section 2.2). The estimates needed to get a Lipschitz continuity prop- 
erty for the composition operator / i— > F o f defined in Sobolev spaces only hold 
for sufficiently regular F's (cf. Section 3, Proposition 2). Runst and Sickel's book 
about Sobolev spaces and Nemystkij operators has been a source of inspiration for 
the treament of the estimates (cf. [32]). The investigation of the Nemytskij operator 
in certain functional spaces has its own interest in functional analysis, independently 
of the special application here (for a study in Lebesgue, Orlicz, Holder or Sobolev 
spaces, see Appell, [12]). Apart from being sufficiently regular, the nonlinearity F 
is chosen to be repulsive to ensure the boundedness of the local solution and by the 
way its global existence. The reason is that the energy of the solution is the sum of 
two terms : one is the part one has even without the nonlinearity (if F = 0) and 
the other one contains the nonlinear part of the energy. The first one is positive 
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(it is directly connected to the norm of the solution in an adapted metric space) 
so the second one has to be positive as well (cf. 2. of the proof of Proposition 1). 
Otherwise, the solution could explode (its norm tends to +00 in a finite time) since 
only the sum of both terms is constant. This restriction leads to what we call a 
"repulsive nonlinearity" which means the dispersion of the energy will increase with 
time. The spatial support of the solution will become larger and larger which keeps 
it from explosion. Segal was one of the first to get interested in the conditions on the 
nonlinearity F which ensure a bounded (and so, global) solution. See [33]. Haraux 
proved that the energy is constant with nonlinearities of the type F(u) = g{u 2 )u 
with some assumptions on g which are equivalent to our conditions on F in the 
second part of Proposition 1 (cf. [19]). The original proof of the existence of global 
solutions of (NLP) with a\ = a 2 > and F(u) = — \u\u\ 2 is due to Jorgens (cf [21]). 

In a second part, the spectral solution of the linearized problem stated in a pre- 
vious work ([10]) is recalled (Theorem 5 of Section 3.2) and causality is proved for 
the linear problem using holomorphy. An adaptation of the Paley- Wiener theorems 
is established which is at the core of the proof (cf. Section 3.2). Then, as in Reed 
and Simon's book (cf. [28]), the Banach fixed point theorem is used to get the 
generalization to the nonlinear problem and due to the assumption F(0) =0 the 
support is respected i.e. the support of F o u is a subset of that of u. 
Note that both parts are quite independent : global existence is not essential for 
causality. 

Our interest in this particular transmission problem dates back to a few years ago. 
In [10], we studied the spectral solution and interpreted the phase gap between the 
original and reflected term as a delay in the reflection. This phenomenon is corre- 
lated to the complementary phenomenon of advanced transmission through a finite 
barrier for low frequency wave packets, theoretically described by Deutch and Low 
in 1993 (cf. [16]). Enders and Nimtz had realized some experiments on superlumi- 
nal barrier traversal before ([18], in 1992). Emig has studied the same situation in 
three dimensions (with waveguides). His conclusion is in accordance with that of 
Deutch and Low : there is an apparent violation of causality. Using the integral 
curves of the energy flow in the space-time he has shown that the maximum of the 
incident wave packet and the maximum of the transmitted pulse are not related by 
this flow : there is a redistribution of the energy contained in the forward tail ([17]). 
Ali Mehmeti and Daikh refined the results of Deutch and Low using the same tech- 
nique of restriction to frequency bands as in our paper of 2003 ([10]). Here again 
causality is not violated by the apparent superluminal effects. At the same time 
Regnier extended this technique to a two-dimensional problem which had been first 
studied from a spectral point of view by Croc and Dermenjian (cf. [14]). 
Transmission problems have been a subject of research for Ali Mehmeti and Nicaise 
since the eighties (cf. [2], [3], [4], [5], [26] and [27]). In 1994 ([6]) Ali Mehmeti gave 
L°°-time decay estimates for the solution of our linearized problem. He found a time 
decay of at least Const.t^ 1 ^ 4 . The reduced decay rate in comparison with the full 
linear case (Const.t^ 1 ^ 2 , cf. Marshall- Strauss- Wainger [24]) seemed to be caused by 
the tunnel effect. But Mihalincic established in 1998 that the decay is Const. t^ 1 ^ 2 , 
being optimal for the time-space points with ^falt > x and Const.t~ m with any 
integer m holds for the other points ([25]). 

In [8], the authors studied the splitting of the energy flow in a star-shaped network 
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for dispersive waves with a same dispersion relation on each branch (a k = a, V k G 
{1 • • -n}). The splitting of the energy flow is independent of the frequency in this 
case. Therefore the phenomenon of delay cannot appear. 

In the future one could study the energy flow for diffraction problems, since Laplace- 
type solution formulae exist (cf. Ali Mehmeti [7] and Rottbrand [30]), as well as for 
scattering problems. 

Concerning nonlinear evolution problems, Haraux's book gives the techniques for 
proving global existence and behaviour of solutions ([19]). Some results about global 
existence of solutions have been established without assuming the nonlinearity to 
be repulsive (cf. [28], [22], [35] and [34]) : dispersion is proved to make up for at- 
tractivity using an L°°-time decay estimate for the solution (cf. [24]). 
As for Kuksin ([23]), he considers a different notion of flow of energy, which is the 
phenomenon of the energy transition to higher frequencies for a nonlinear wave equa- 
tion under n-dimensional periodic boundary conditions in a Sobolev phase space. 
Note that the nonlinearity keeps the energy from staying in a frequency band. 



2 Existence and uniqueness of solutions 

2.1 Abstract setting 

Let us consider the data as in the introduction and define the linear operator 
(A, D(A)) as follows : 

D(A) = {ve nLi H 2 (n k )/v satisfies (T ), (71)} 

A: v E D(A) ^ Av := f-c 2 ^| + a k v k ) in H := njU L 2 (n k ) 

\ dX /fce{l;2} 

and given u(t) = (ui(t, •), u 2 (t : •)) and F(u(t)) = (F(ui(t, •)), F(u 2 (t, •))), the ini- 
tial boundary value problem (NLP) (see the introduction) can be rewritten as the 
following problem denoted by (P) : 



(P){ 



' u e C 2 (IR, H),u{t) e D(A), \/teIR + 

d 2 u 

(t) + Au(t) = F(u(t)) in H,VteI& 



dt 2 

«(o) = ^(0) = o 



where $ = (f u / 2 ) G D(A) with f k :tt k ^ IR for k G {1; 2}. 

Note that for the sake of simplicity, the initial velocity Ut(0) is chosen to vanish but 
putting Ut(0) = would not change the problem fundamentally. 
The spectral theory of the linearized problem (F = 0) is studied in [7], in Theo- 
rem 1.5.1. Here dj = 1, V j G {1 • • -n} and for us the c/s and a/s have to be 
interchanged. The operator A is proved to be self-adjoint and the existence and 
uniqueness of a solution of the linearized problem (LP) ((LP) is (P) with F = 0) 
in C 2 (IR + , nLi L 2 (n k )) n C\IR + , D(VA)) n C°(IR + , D{A)) is stated in [7] (in the 
context of the analysis of the asymptotic time behaviour of the solution). Recall 
that D(VA) := {v G Ul=i H\n k ) /v satisfies (T )}. 
As Reed and Simon do in [28], one can rewrite (P) as (Q) : 

(Q){ ^) +^(0 = <%(*)) in H 
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where tp(t)= ( du j and J(<p(t))= ^ F (^)) ) and A = * ( _ ^ q ) 

The operator .4. is proved to be self-adjoint and closed on its domain of definition 
D(A) ® D(VA) (cf. the proposition of Section X.13 in [28]). 



2.2 Reed-Simon's Theorems for the existence of solutions 

Let us recall the three main theorems given in [28]. Reed and Simon's proofs are 
all based on the contraction mapping principle applied to the mapping contained in 
the reformulation of Problem (P) as the integral equation 

<p(t) = e -^Vo+ f e-^- s)A J{ip{s))ds 
Jo 

For any <p G D(A), the norm of ip — (u,v) is defined by Reed and Simon as 

y\\ = (\\VAu\\i + \\v\\i)^. 

Thus, in our case, A<p, J(<f) and AJ((p) are defined as 

\\A<p\\ = (\\VAv\\l + \\Au\\iy/2 
\\J( V )\\ = \\Fou\\ 2 
\\AJ(y>)\\ = \\VA(Fo U )\\ 2 

Note that, in our case, u(t) has two components so ||ii|| 2 has to be understood as 

/ \u k {t,x)\ 2 dx 
\k=i jQk 

In the following three subsections, we recall general results on the existence and 
smoothness of the solutions of (Q), where A is any self-adjoint operator and J any 
mapping from D(A) to D(A) (see in [28] for the proofs). 



2.2.1 Local existence 
Theorem 1 

Let A be a self-adjoint operator on a Hilbert space H and J a mapping from D(A) 
to D(A) which satisfies J(0) = and : 



(h l ) \\j(<p)-jm<c(\\<p\\,w\\)-\\<p-i>\\ 

(iff) \\A(J(<p) - < C(M, \\Ap\\, Ul ll^ll) • \\Ap - A1>\\ 

for all ((f,ip) G D(A) 2 where each constant C is a monotone increasing (everywhere 
finite) function of the norms indicated. Then, for each ip G D(A), there is aT > 
so that Problem (Q) of Section 2.1 has a unique continuously differentiable solution 
for t G [0; T). For each set of the form {ip / \\(p\\ < a, \\A(fi\\ < b}, T can be chosen 
uniformly for all ip in the set. 



Remark 1 

Note that the assumptions (H ) and {Hi) of Reed and Simon (see [28]) are implied 
by J(0) = ; (Hq) and (H^). Moreover (H ) and {Hi) make the condition J(0) = 
necessary. Hence the simplifications we have made here. 
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2.2.2 Local smoothness 



Theorem 2 

(a) Let A be a self-adjoint operator on a Hilbert space H and n a positive integer. 
Let J be a mapping which takes D(A*) to D(A>) for all 1 < j < n (for allO < j < n) 
and which satisfies J(0) = and : 



{Hf) \\A{m - < cm, mvii, i wid • \\a?<p - a^\\ 



for all ((p,ip) G D(Ai) 2 where each constant C is a monotone increasing (every- 
where finite) function of all its arguments. Then, for each tpo G D(A n ), n > 1, 
there is a T n so that Problem (Q) of Section 2.1 has a unique solution ip(t) for 
t G [0,T n ) with ip(t) G D(A n ) for all t G [0,T n ). For each set of the form 
{ip j ||«4V|| < %; j = 0, n}, T n can be chosen uniformly for all ipo in the set. 

(b) In addition to the hypotheses in (a), assume that for each j < n, J has the 
following property : if a solution ip is j times strongly continuously differentiable 
with y2 (fe) (t) G D(A a ~ k ) and A n ~ k (p ( - k \t) is continuous for all k < j, then J((p(t)) 
is j times differentiable, d j J(ip(t))/dt j G DiA 71 ^' 1 ), and A n ~ j ~ 1 d j J(ip(t))/ dt j is 
continuous. Then the solution given in part (a) is n times strongly differentiable in 
t and d j J((p(t))/dt j G D(A n - 3 - r ). 

2.2.3 Global existence 
Theorem 3 

Let A be a self-adjoint operator on a Hilbert space H and n a positive integer. Let 
J be a mapping which takes D{A^) to D(A') for all 1 < j < n (for all < j < n) 
and which satisfies J(0) = and : 



(H>) \\A*J{<p)\\ < C(|MI, II^-VII) • MVII, 1 < j < n 

{Hf) \\J?(J(<p) - J(M\ < C(hl HVII, MVII, ■ \\Ai<p - Aty\\ 



for all {<p,ip) G D(Ai) 2 where each constant C is a monotone increasing (everywhere 
finite) function of all its arguments. Let tp G D(A n ) and suppose that, on any finite 
interval of existence, the solution ip(t) guaranteed by part (a) of Theorem 2 has the 
property that \\(p(t)\\ is bounded from above. Then there is a strongly differentiable 
D(A n ) -valued function (pit) on [0; +oo) that satisfies 



Further, if J satisfies the hypotheses of part (b) of Theorem 2, then <p(t) is n-times 
strongly differentiable and d^ip(t)/dP G D{A n '^). 



2.3 Required estimates 

In this section, we will apply the last three theorems to our particular operator A and 
mapping J that is we will first interpret the required estimates of Theorem 3 in terms 



< j < n 
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of a Lipschitz-type continuity for the nonlinearity F (Proposition 1) and establish 
the required regularity properties for F to satisfy those estimates (Proposition 2). 
At last we will deduce the assumptions on F necessary for the existence of a global 
solution to Problem (P) (Corollary 1). Analogous results can be found in Haraux's 
book (A. II. Lecture 3 of [19]) : in particular he proves that the energy is constant for 
nonlinear wave equations with a nonlinearity of the type F(u) = g(u 2 )u under some 
assumptions on g in dimension N (g' has to lie in Lf£ c (IR + ) in dimension one). Note 
that we restrict ourselves to the proof of the one-time strong differentiability of the 
solution (n — 1 in Theorem 3). To get a higher regularity, nonlinear compatibility 
conditions have to be added to the transmission conditions (T ) and (Ti) (cf. [4]). 

Proposition 1 

1. For the conditions J(0) = 0, (Hq) and (H[) to be satisfied, it is 

• necessary that F(0) = 

• sufficient that 

F(0) = 

||F O u - F O y\\ 2 < C(||u|| 2 , ||Ux||2, \\v\\2, \\vxh) ■ \\VAu - VAv\\ 2 

\\VA(F ou-Fov)\\ 2 < C(\\u\\ 2 , \\u x \\ 2 , \\v\\ 2 , \\v x \\ 2 ) ■ \\A(u - v)\\ 2 

for all (u,v) G D(A) 2 and where C is a monotone increasing function of 
all its arguments. 

2. Assume F e C°(IR) and suppose that for any (u,v) G (IljLi L 2 (Qk)) 2 , there 
exists C(||«||, H^H), monotone increasing function of all its arguments such that 

\\Fo U -Fov\\ 2 < C(|M| 2 , ||u|| 2 ) • ||«-u|| 2 

IfG, defined on IR by G(w) = f$ F(£)d£, takes only non positive values and if 
it is such that Gou belongs to L x (f2) for u G H 2 (Q) (fl being an interval of IR, 
possibly unbounded), then the unique local solution of Problem (P), guaranteed 
by the first part of the proposition in combination with Theorem 1 and denoted 
by ip, has the property that ||y?(£)|| is bounded from above on any finite interval 
of existence. Thus, due to Theorem 3, it defines a global solution of (P) i.e. a 
strongly differentiable D {A) -valued function on [0; +oo) that satisfies Problem 
(P)- 

Proof. 

Let us recall that the norm of tp — (u,v) is defined as \\tp\\ = (\\VAu\\l + |M|!) 1/2 
for any tp G D(A). And those of Atp, J{<p) and AJ(tp) are defined as 

11^11 = (\\VAv\\l + \\Au\\iy/2 
\\J( V )\\ = \\Fou\\ 2 
\\AJ(tp)\\ = \\VA(Fo U )\\ 2 

1. • The necessary condition is clearly required : J(0) = and J{tp{t))— 
( FW)) ) imP ' y Fi0) = °" 
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Now, if is still denned as (p = (u, v) and ip as ip = (u, v). So, if 

||F o u - F o u\\ 2 < C(||u|| 2 , \\u x \\ 2 , \\u\\ 2 , Wuxh) ■ \\VAu - VAu\\ 2 
then \\F o u — F o u\\ 2 

< C(\\u\\ 2 , \\u x \\ 2 , \\u\\ 2 , \\u x h) ■ (\\VAu - ^~Au\\ 2 2 + ||u - v\\ 2 2 ) 1/2 
Now, the norm of \fAu is defined by : 

2 

1 1 V^4i* 1 1 2 = (VAu; y/Au) = (Au,u) = \\ u k\\ 2 + a fc ||w fcrr || 2 

fc=i 

so \WxW2 < max(l/ai ; 1/02) || 2 < max(l/ai; l/a 2 )IMI an d since C 
is a monotone increasing function of all its arguments, C(||m|| 2 , Hw^lh) is 
lower than C"(||<^||), with C"(||<^||) a monotone increasing function that 
will be denoted by C as well. 
Likewise C(||u|| 2 , IKH2, INI2, \\uxh) < C(IM|, 



fr- 



aud (\\VAu - VAu\\l + \\v - d\\j) 1/2 
So the first condition given in the first part of the theorem, has been 
proved to be sufficient for (Hq) to hold. Idem for iff. 

2. If u is the solution of Problem (P), the energy is defined as 

E ^) = \T. I (\(VAu(t,x)) k 2 + \u k , t (t,x)\ 2 -2 (Gou k )(t,x))dx 



x)dx. 



that is: E{t) = \\\<p{t)f-Y, I (Gou k )(t, 

Note that Gou k belongs to L 1 (Q fe ) by assumption (u is the solution so u k (t, •) 
lies in H 2 (Q k )). Thus, since G takes only non positive values, 

^ll^)H 2 <^ll^)l| 2 -E / (Gou k )(t,x)dx 



i.e. -y(t)f<E(t). 

Now let us show that the energy is independent of t to get 7:||y?(£)|| 2 < E(0) 
i.e. a bounded norm of ip(t) on its interval of existence. 

The way it will be done is the same one as in Haraux ([19], p. 23). His 
nonlinearity is of the form F(u) = u-g(u 2 ) and his conditions on g in dimension 
one are equivalent to ours. 

So first of all, the first two terms in E(t) are clearly differentiable. The reason 
is the same one as in Reed and Simon's book : if is strongly differentiable as 
a D(A)-valued function. It means in particular that : 



A [ u(t + h)-u(t) _ ut{t) " 



h 







' u t {t + h) - u t {t) 
h 



- u tt (t) 
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as h — > 0. 

There remains to deal with the third term containing the nonlinearity : 
Since F G C°(IR) and since G is defined on IR by G{w) = J™ F(£)d£, G 
is differentiable and G' = F. Now, for any (w,l) G IR 2 , there exists c in 
[min(u>; w + I); max(u>; w + I)] such that, 

G(w + I) - G(w) = lF(c) 

Applied to w — u k (t, x) and I — u k (t + h, x) — u k (t, x), it reads : 

G(u k (t + h,x))- G(u k (t, x)) = F(y%(t, x)) ■ (u k (t + h,x)- u k (t, x)) 

with y k (t, x) G [mm(u k (t, x); u k (t + h, x)); max(u k (t, x); u k (t + h, x))]. 
Note that we should write F(y h (t,x)) k with y h (t,x) = (yi(t,x);y%(t,x)) since 
F has been defined on the product space H . We will always keep this notation 
in the following. Now 



; G(u k (t +h , x)) - G (u k (t, x) ) _ F{uk{t)X))ukt ^ x) 



dx < 



Ia k ^{[G(u k (t + h,x)) -G(u k (t,x))] - F(u k (t,x))(u k (t + h,x) -u k (t,x))} 



dx 



+ J Qk \F(u k (t, x)) ( "»(«+M)-u»ft») _ UkAtj x) j 



dx 



Then, due to Holder's inequality 



E 2^ | G { u k{ t +h , x)) -G(u k{ t, x)) _ x))uk t ^ x) 

(H k \\F o y h k (t) - F o u k {t)\\ 2 ) 



< 
< 



u(t+h)-u(t) 
h 



dx 



+ \\F(u(t))h 



(E*c(||^(0l|2;IM*)ll2)-||^(*)-« fc (0ll2) 



+\\F(u(t))h 



u(t+h)-u(t) 
h 



u(t+h)-u(t) 
h 



u(t+h)-u(t) 
h 



- Ut(t) 



~ U t (t) 



Now "(*+ ft )-"( f ) tends to u t (t) in L 2 since u G C\IR,H) (it is even of class 
C 2 ). So there remains to prove that \\Foy^(t) —Fou k (t)\\2 tends to zero with 
h to conclude that the latter estimate tends to 0, which means that the third 
term is differentiable and the derivative of E(t) is then 

2 r 

E '(t) = E / u kAt, x) (u k , tt (t, x) + (Au) k (t, x) - [F(u(t, x))] k ) = 

due to Equation (E k ). 
First of all, it holds classically : 

min(/; I') 

< 

max(7; I') 



= hi + i'-\i-r\) 
= -(i + i' + \i-i'\) 



9 



Then 

\yk(t,x)\ 2 < ^{u{t,x) +u{t + h,x) + \u(t,x) -u(t + h,x)\) 2 
< A(\u(t,x)\ 2 + \u(t + h,x)\ 2 ) 

and y%(t, ■) belongs to L 2 (VL k ) since both u{t,-) and u(t + h, ■) are in L 2 . 
Moreover ||y^)|||<4(|Kt)||I +11^ + ^)111). 

Now y k (t, x) is superior to m := ~(u(t, x) + u(t + h, x) — \u(t, x) — u(t + h,x)\) 

1 

and inferior to M := —(u(t, x) + u(t + h, x) + \u(t, x) — u(t + h, x)\). 

Then \y%(t, x) — u(t,x)\ is inferior to the maximum of m and M. So it is 
bounded by (u(t + h, x) -u(t,x)). Thus \\y%(t) -«(*)||2 < \\u(t + h) - u(t)\\ 2 
and, since u G C°(IR, H), u(t + h) tends to u(t) in L 2 . 
Then lim^o \\y£(t) - u(t)\\ 2 = 0. ■ 

Conditions on the nonlinearity of the problem i.e. on the function F can be deduced 
from the latter proposition. Runst and Sickel's book (cf. [32]) has been a source of 
inspiration for the following and especially for the next two proofs. Using the Fatou 
property of the Sobolev space H 2 ((0; +oo)) is their key idea (cf. [32], p. 15). 

Notation 1 

The composition operator T F : j ^ p Q j satisfies Property Cip(B, B') if there 

exists a constant C > 0, monotone increasing function 0/H/H2, IMI2, H/xlh, Hfl^lh; 
such that for any (f,g) in B 2 C (L 2 ) 2 

\\Fof-Fog\\ BI <C(\\f\\ 2 ,\\gh,\\f x \\ 2 ,\\g x \\ 2 )-\\f-g\\ B 

where || ■ || 2 is the norm in L 2 . 

Proposition 2 

Let F : IR — > IR be a continuous function satisfying F(0) = 0. Denote by Tp the 
composition operator Tp '■ f 1— > F o f 

Then, for Tp to satisfy Property Cipi^H 1 ((0; +00)), L 2 ((0; +00))) (resp. Cip(H 2 , H 1 )) , 
it is : 

• necessary that F e H 2 ' loc ((0; +00)) 
(respectively F e H^ loc ((0; +00))) 

• sufficient that F G C 2 {IR, IR) (respectively F G C l {IR, IR) ) 
Proof. 

Runst and Sickel's proof is adapted since the estimate we need is a weaker result as 
theirs : it contains a constant depending on four parameters, so we will explicit this 
function and check its monotony. 
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1. Necessary condition : Choosing a sequence of functions fu in S(IR; IR) such 
that I— > /m(^) = £ when |x| < M leads to T F (f M )(x) — (F o /m)(^) = 
F(x) when |x| < M i.e. T F (f M ) = F\^ M . M y Thus for T F (f M ) to belong to 
# 2 ((0; +00)) (respectively ^((O; +00))), F\[_ M . M] has to lie in H 2 ((0; +00)) 
(resp. ^((O; +00))) for any M in IR i.e. F has to lie in H 2 > loc ((0; +00)) (resp. 
//^((Oj+oo))). 

2. Sufficient condition : the aim is to prove that, if F G C 2 (IR), V(/, (7) in 
(/J 2 ((0;+oo))) 2 ) 3C(||/|| ) ||^| ) ||M| ) ||^||)>0 

llFoz-Fo^i^^cdi/iUbiuiMUi^iDH/-^^ 

i.e. Tp satisfies Property Cip(H 2 , H 1 ). Since the arguments for CiplH 1 , L 2 ) 
are strictly analogous, this case is not treated separately. 

• First step. 

Let us first prove that Property £ip(C™(IR)C)H 2 ((0; +00)), H 1 ((0; +00))) 
is satisfied for F G C 2 (IR). 

The norm ||F o / — F o g\\m(n) is equivalent to ||F o / — F o qWl'^q) + 
\\d x (Fof — F o g)\\ L 2(Q} for any domain Q. So the aim is to find estimates 
for \\F0f-F0 g\\ L2{n) and ||(F o /)(/') - (F o g)(g')\\ L 2 (n) assuming 
that (f,g) G (C c °°(/i?) n H 2 ((0; +00))) 2 and F G C 2 (/#). 



(a) \\F o f - F o g\\h m+oo)) = J +o ° \(F o /)(*) - (F o ^)(x)| 2 ^ and 
|(F o f)(x) - (F o ^)(x)| = \F(f(x)) - F(^(x))|. Thus, for any 
x G (0; +00) : 



\{{F o f) - {F o g)){x)\ < 



sup I if 

\y\<max{\\f\\ L ca;\\g\\ L oo) 



\f(x)-g{x)\ 



So ||F o / - F o ^|| i2((0;+ oo)) < M(F, /, • 11/ - <?|| L2 ( ( o ;+ oo)) with 
M(FJ,g)= sup \F'(y)\ 

\y\<max(\\f\\ L oo;\\g\\ L oc) 

(b) ||(F'o/)./'-(F'o y ).^||i 2 

< 2||(F' o /) ■ (/' - g')\\ 2 L2 + 2||((F o /) - (F' o g)) ■ ^||| 2 

<2/o +00 (^ , (/W)) 2 (/ / -^) 2 (^)^ 

+2 / +oo (F(/(:r)) - F'(g(x))) 2 (g') 2 (x)dx 

Now the first integral is bounded from above by 



sup \F' 

■\<WfU°° 



And the second one by 



sup 

Jj/l^maxdl/HioojIlffHioo) 



\F"{y)f 



Wf-g'Wh 



\\{9)'\\h-\\f-g\\ 



Thus, || (F' of)-f- (F> o g) ■ g'\\ L2 < M'(F, f, g).\\f - g\\ mm+oo)) 
with 

M'(F, f,g) = Const- max ( sup M < i|/!|LOO |F'(y)|; ||^|| i2 . 

SU P|l/|<max(||/|| LO c;||p|| LO c) \ F " {v)\ ) 
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In this first step of the proof, we have then proved that, if F E C 2 (IR), 
\\Fof-Fo #||ffi((0;+oo)) < ^ ) (^/ ) 5 , )ll/-5 , l|// 1 ((0;+oo)), for any (f,g) in 
(C™(IR) n H 2 ((0; +oo))) 2 . We have denoted by D, the function of F, f 
and g : D(F, f, g) = Const ■ (M(F, f, g) + M'(F, f, g)) where M(F, f, g) 
and M'(F, f,g) have been defined previously. 

Second step. 

Using density and the Fatou property of the Hilbert space if 1 ((0; +00)), 
we will deduce from the first step that Property Lip is also satisfied for 
any (/, g) E (# 2 ((0; +00))) 2 and F E C 2 {IR, IR). 

Since C™(IR) n H 2 ((0; +00)) is dense into H 2 ((0; +00)), there exists a 
sequel (f k ) included in C%°(IR) fl H 2 ((0; +00)) such that f k converges to 
/ in H 2 . Yet L°° is continuously embedded in H 1 so f k converges to / in 
L°° and so, almost everywhere as well. The continuity of F implies that 
F o f k tends to F o / a.e. as — ► +00. 

Now, for any k E IN, F and f k are continuous and, for any bounded inter- 
val Q of IR, the restriction of / to Q admits a continuous representative 
as an H 2 function. So F o f k and F o f belong to L 1,/oc ((0; +00)) and, 
due to Lebesgue's theorem, F o f k converges weakly to F o f in S' as 
k — >■ +00. 

The same holds for g which can be considered as the limit of a sequel g k . 
By assumption, for any k E IN, F E C 2 (IR) and (f k ,g k ) belonging to 
(C c °°(ffi) n H 2 ((0; +00))) 2 , then there exists a constant D(F, f k ,g k ) > 
such that ||Fo/ fc -Fo0 fc || H i < D(F,f k ,g k ) ■ \\f k -g k \\m for any k E IN. 
Since f k (resp. g k ) converges to / (resp. g) in H 2 , Mk k -^ +0 o\\fk ~ 9k\\w 
is equal to ||/ — g\\m- Now 

D(F, f k , g k ) = Const ■ (M(F, f k , g k ) + M'(F, f k , g k )) 

where 

' M{F,f k ,g k ) = snp lyl ^ maxmUoc . llgkUoo) \F'(y)\ 
M\F,f k ,g k ) = v / 2-max(sup M < )|/fc||i0o \F'(y)\; \\(g k )'\\v 

[^Pll/l^moxdl/fclUoojHsfcllioo) \ F "(V)\ ) 

Since f k converges to / in L°°, \\f k \\L°° tends to ||/||l°° and so, M(F, f k , g k ) 
tends to M(F, f, g) as k tends to 00. On the other hand, g k converges to g 

in H 2 implies that (g k )' converges to g' in L 2 . So lim fc i+OQ M'(F, f k , g k ) = 

M'(F,f,g). 

Thus hm fe _ +00 ||F of k -Fo 9k \\ Hl < D(F, f,g)-\\f- 9 \\h*. 
Now, the Hilbert space i? 1 ((0; +00)) is such that S C H 1 C S' and 
has the Fatou property (cf. [32], p. 15, for the definition of the Fatou 
property). Thus, V / E H 2 , (F o / - F o g) E H 1 and 

\\Fof-F og\\ m < Const ■ D(F, f,g)-\\f- g\\m 

It can also be written : ||F o / - F o g\\ H i < £?'(||/|| 2 , ||^|| 2 , ||/'|| 2 , ||gr'|| 2 ) • 
11/ - g\\m with 

L>'(F,||/|| 2 ,y| 2 ,||/'|| 2 ,|| 5 '|| 2 )=C^-( \F'(y)\ 
+ max(sup| 2/ |< Cst .|| / | |wl \F'(y)\; \\g'\\ 2 ) ■ [sup^!^.^^^^^^^ \F"{y)\] ) 
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since H 1 and H 2 are continuously embedded in L°°. 
To finish with, there lacks to make sure that D'(F,f,g) is a monotone 
increasing function of its arguments : by definition, H/H/fi is a monotone 
increasing function of ||/|| 2 and ||/'||2- Moreover M h- > sup^ <M \F'(y)\ is 
a monotone increasing function so is D'. u 

Remark 2 

In the proof, we have also stated that, if F belongs to C 1 (IR), then for any (u,v) 
in (nik=i L 2 (Q k )) 2 , there exists C{\\u\\2, \\v \\ 2 ), monotone increasing function of its 
arguments such that 

\\F ou — Foi;||2<C(||«||2, || ^ || 2) " || M — f H2 

which is needed in the second part of Proposition 1 to get the boundedness of the 
solution and so, its global existence. 

Corollary 1 

Let F G C 2 (IR,IR) such that F(0) = 0. // G : w i-> fif F(£)d£ takes only non 
positive values on IR and if G o u belongs to for u G L 2 (Q) (Q being an 

interval of IR, possibly unbounded) , then the unique local solution of Problem (P) 
guaranteed by the first part of Proposition 1 defines a global solution of (P) i.e. a 
strongly differentiable D (A) -valued function on [0; +00) that satisfies Problem (P). 

Proof. 

Due to Proposition 1, we only have to check that the following sufficient conditions 
are satisfied : 

f ||P o u - F o u|| 2 < C(|H| 2 , |K|| 2 , 1Mb, IKH2) • \\VAu - VAv\\ 2 
\ \\VA(Fou-Fov)\\ 2 < C(|H| 2 ,||u x || 2 ,||u||2,K||2) • \\A(u-v)\\ 2 

for all (u,v) G D(A) 2 and where C is a monotone increasing function of all its 
arguments. 

Since the sufficient condition of Proposition 2 is satisfied, it holds for all (w, v) in 
D{A) 2 : 

f \\Fou - F ov\\ 2 < C(||u|| 2 , ||u x || 2 , \\v\\ 2 , \\v x \\ 2 ) ■ ||u-u||ffi 
\ \\F o u - F o v\\m < C(\\u\\ 2 , \\u x \\ 2 , \\v\\ 2 , \\v x \\ 2 ) ■ \\u - v\\ H * 

So the first estimate reads : 

ELl \\F OU k - F O V k \\ 2 2 < C 2 ELl ||«fc - v k\\l + Efc=l ll M M ~ V k,x\\l 



<C 2 -M- 
<C 2 -M- 



^2k=l\\ U k v k\\ 2 ^2k=l a k\\ u k,x v k,x\\ 2 

\\fAu - \fAv\\\ 



i.e. ||Pow-Pod| 2 < C • \[M ■ \\VAu\\ 2 with M = max (— ; — ; lV Since M does 

Vai a 2 J 

not depend on u nor on v, C ■ \f~M is still a monotone increasing function of ||m|| 2 , 
IK|| 2 , |H| 2 and \v x \ 2 . 

Moreover the second estimate given by Proposition 2 reads : 
ELi \\Fou k -Fov k \\l + \\(Fou k ) x - (Fov k ) x \\ 2 
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< c 2 

Then 



Efc=l \\ u k v k || 2 ~l~ Efc=l \\ u k,x ^k,x || 2 Efc=l II u k,xx Vk y xx\\ 2 



\\y/A(Fo u- Fov)\\l 

= \\(Fou 1 -Fov 1 )\\l + a 1 -\\(Fo Ul -Fov 1 ) x \\l + \\(Fou 2 -Fov 2 )\\ 2 2 

+a 2 ■ \\(F o u 2 - Fo v 2 ) x \\l 
< a 2 • (ELi \\F ou k - F o v k \\l + || (F o u k )x - (Fo v k ) x \\fj 

<C 2 -a 2 -M- [ELi IK - V k\\l + Efe=l \\ U k,xx - Vk,xx\\l 

+C 2 -a 2 - M -j:t =1 \\uk-v k \\ 2 
<C 2 -a 2 -M- \\Au\\ 2 + C 2 ■ a 2 ■ M x 2 (sup fcg{1;2} K 2 ) ■ \\Au - Av\\ 2 

where K k is defined through an Ehrling-Nirenberg-Gagliardo (interpolation) inequal- 
ity (cf. Adams [1], Th 4.14). With m = p = 2 and j — 1, it reads for a fixed e > 
and any < e 

II II 2 -C K II II 2 j_ K — 1|| II 2 
|| w fc,a;||2 — -^A; e fc|| M fc||2 ' ^k^k y^k,xx\\ 2 

Choosing = e = 1 leads to the latter estimate. We have then proved that 

\\[A(F ou-F ov)f 2 < Const - \\Au - Av\\ 2 

And the constant is a monotone increasing function of all its arguments : ||w||2, 
\\uxh, 1Mb and ||f x ||2- ■ 

Example. If the nonlinearity F is chosen to be defined by F(u) = —Aw 3 with 
A > 0, it satisfies all the conditions required in Corollary 1 : F(0) = is clear as 
well as the regularity (F G C 2 (IR)). 

4 

W 

G : w i— > Jq F(^)d^ takes only non positive values on IR : in fact, G(w) = — A— 

and w 4 > for any w G IR. At last, since G is continuous and since any u G H 2 (Q) 
is also continuous on Q, G on lies in L 1,loc (Q). As for the behaviour at infinity (if Q 
is unbounded), a continuous L 2 function is in L 4 . 



3 Causality of the support of the solution 

3.1 Spectral solution formula of the linearized problem (LP) 

The aim of this section is to rewrite the spectral solution formula of the linearized 
problem already studied in [9] and [10], assuming that the initial condition is com- 
pactly supported in the first branch (0; +oo). Let us recall that the expression for 
the unique solution of the problem is given by the following theorem : 

Theorem 4 

Assume that f k '■ Q k — ¥ IR with k G {1;2} ; initial data of Problem (LP), described 
in the introduction, are such that (f\, f'2) G D(A). 

Then the restriction to Q\ of the unique solution of Problem (LP) (such that (ui,u 2 ) 
belongs to C 2 (IR + , IlLi L 2 (tt k )) n C\IR + , D(y/A)) n C°(IR + , D(A)) ) is given, for 
(t,x) in IR + * xfli, by : 
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+ ^ Wo) cosM)S ( ^j^W) (/ +O ° /i(«)e-^^)^)d«)) 

+^ Wo) cosM)3 ( ^g^^ (/ +O ° h(u)e-^-«^du)) <L> 



R, D(A), D(y A) are introduced in Section 1. 
Recall that, for real u and j G {1; 2} 



' Oi — w 2 



— — if u; 2 < a, 
cf 

^2 lf W 



Note that the solution is the sum of three terms : the original term, the reflected 
one which appears due to the discontinuity in the potential and the last one, called 
"transmitted term" , which contains the information on the transmission of the signal 
from the second branch to the first one. Since our point is to study causality, we 
will restrict ourselves to a compactly supported initial condition with support in 
the first branch and the third term will vanish. Rewriting the solution in terms of 
Fourier transforms, we get : 

Theorem 5 

Assume that fj, : Qj, — > IR with k G {1; 2} ; initial data of Problem (LP), described in 
the introduction, are such that / 2 ) G D(A), f\ is compactly supported in (0; +oo) 
and fi = 0. 

Then the restriction to Qi of the unique solution of Problem (LP) (such that (ui,u 2 ) 
belongs to C 2 (IR + , n*=i L 2 (n k )) n D(^fA)) n C°(IR + , D(A))) is given, for 

(t,x) in IR + * x n 1 , by : 



ui(t,x) = -F^Lx cos(v/a^ + c 2 uj 2 t)J r f 1 (uj) 

7T 

1 



71 



■-1 
a) i— >x 



I I I £~S4.\ I uc-V c*u* - a 2 + a t 

cos (v o-i + c2ujH ) — i / 2 2 =r= 

\uc + v c ^ — a 2 + a 1 



2, ,2 



where the complex square root has been defined such that ^c 2 uj 2 — (a 2 — ai) = 



a2— ai 
„2 



c 2 u 2 — (a 2 — ai) if uj > 
Wa 2 - ai - c 2 u 2 if -J^^<co< 



— yc 2 uj 2 — (a 2 — ai) if ^ < — 
Likewise, for (t,x) in IR + * x Q 2 , u 2 (t,x) = 

i ^ i ^/ c 2 w 2 — a2+ai^ a; / 



«2 — Ql 
„2 



1 Z" 00 / i 

/ cos ( J ai + c 2 uj 2 t)e 

2% J-<x 



9.1, ir 



k ujc + \Jc 2 uj 2 — a 2 + ai 
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The functions of u denoted by Cr and T 1,2 and defined by 



C r (uj) :-- 



ujc — \Jc 2 uj 2 — a 2 + a[ 
ujc + v 'c 2 ui 2 — a 2 + ai 



and T ' (u;) : 



2wc 



cjc + \J c 2 ui 2 — a 2 + ai 



are called reflection coefficient and transmission coefficient (from the first branch to 
the second one) respectively. 



Proof. 

First of all, in the integral defining u\, uj is changed into u/ = £i(c<; 2 
i.e. = V^i + c 2 uj' 2 and do/ = c 2^ 2 ) = ^Mw)' Tlie domain of integration 
((V^ii +°°)) thus becomes (0; +00) with the reflection coefficient taking different 
values depending on the belonging of uj to 0; 
It holds : 



a,2 — a\ 



or 



(12 — ai ■ 



«i(*,a;) 



1 



3? 



7T 



cos( v /a 1 + c 2 cu 2 t)e-^/;(cj)^ 



7T 

+ -» 

7T 



cos(\/ai + c 2 uj 2 t)e l 



lUOC 



a 2 — a x 



c 2 u 2 



iujc 



+ \fa 2 — a x — 



C 2 UJ 2 



r+oc 



cos( Jai + c 2 oo 2 t)e z 



'ujc — \Jc 2 uj 2 — a 2 + ai 
ujc + \Jc 2 uj 2 — a 2 + ai 



j J r f 1 {-uj)du 
Tji (—u)duj 



In the second step of the proof, is replaced by (z + z)/2 and is changed into 
— uj in the conjugate term. The complex square root of the reflection coefficient is 
chosen to be analytic in the upper half-plane as in Deutch and Low's paper (cf. [16] 
and the above definition of the complex square root 5). 

As for u 2 , its expression for a compact supported / is reduced to one term analogous 
to those in U\ (cf. Theorem 5) : u 2 (t,x) = 

Woo) ( ^fZ^K^ {Jo - f{u)c^>-^du)) duj 

Then the same type of modifications are applied to get the formula of Theorem 
5. ■ 



3.2 A Paley/ Wiener- type theorem for a function holomor- 
phic in the upper half-plane 

The aim is to generalize a classical Paley- Wiener Theorem which gives explicit 
conditions on L 2 (/i?)-functions to be the inverse Fourier transforms of compactly 
supported functions. This theorem, cited and proved in Rudin ([31]), requires ana- 
lyticity on the whole plane as well as an exponential estimate and a L 2 -behaviour 
on the real line. 

Our solution formula for the linearized problem given in Section 3.1 involves the 
inverse Fourier transforms of functions which satisfy those requirements except the 
analyticity on the whole plane. The square roots appearing in both the reflection 
and transmission coefficients keep them from defining entire functions. They can 
only be analytic in the upper half-plane (for example) 11+ = {z e (D/Q^z) > 0}. 
And we will show that the difference with an entire function is that the support of 
the inverse Fourier transform of such a function is not compact but only bounded 
from above. 
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So our goal here is to adapt the proof of Theorem 19.3 of [31] which is based on 
Cauchy's Theorem. 

First of all, we need another representation of the square root \/c 2 u 2 — a 2 + ai 
to avoid the branch cuts to lie on the real line (it will be useful to apply Cauchy's 
Theorem in the proof). Writing, as Ali Mehmeti does in [7] 

a \fz := yj~2|e^ arga(z) where z = re iarga(2) , r > 0, &rg a (z) G [a; a + 2n) 

we choose y/c 2 uj 2 — a 2 + a\ to be ^ n / 2 y/cu> — kj • (~n/2V CUJ + kj with k = ^a 2 — a[. 
Thus the branch cuts become 7(±A;) = {±k — iy/y > 0}. This representation does 
not change the value of the square root on the real line. In particular it is still 
continuous on the closure of H+ as constructed in the last subsection. 

Theorem 6 

Let A and C be two positive constants and g : € — ► € be such that : 

• g is analytic everywhere except for the two cuts : 
7(±ib) = {±k - iy/y > 0} where k G IR + * 

• \g(z)\ < Ce A \ z \ for any z G 11+ = {z G C/5S(z) > 0} 

• IlR \9(x)\ 2 dx < +oo 

Define the function G on IR by G = T g where T is the L 2 -Fourier transform. 
Then the support of G is a subset of (—A; +oo) and G G L 2 (—A] +oo). 

The following proof is an adaptation of Rudin's one (cf. [31]). 

Proof. 

For e > and any real x, let g e be defined by g t = g(x)e~ e ^. The aim is to show 
that, for any real t < —A : 

\imj iR g e (x)e- u *dx = (*) 

Since \\g — gA\ 2 — > as e — > 0, Plancherel's Theorem implies that the Fourier 
transform of g e (denoted by G e ) converges towards the Fourier transform G of the 
restriction of g to the real axis. Then a subsequence G en converges almost every- 
where towards G and, since lim n _ +00 G en vanishes outside [— A; +oo) (due to (*)), 
so does G. 

So proving (*) for any real t < — A is sufficient. 

For any real a G [0; 7r], we define like Rudin : 

r Q (s) = se ia , < s < +oo and Tl a = {w/$t(we ia ) > A} 
And for w G Tl a : 

$ Q H = / g(z)e-™dz 
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It is clear that <3> Q is defined on IT a (g is continuous on the closure of n + so it is on 
T a and \g(z)e~ uz \ < Ce' 3 ^^-^). It is analytic on Tl a for a e (0; tt) (cf. [31]). 
Moreover it holds for any real t : 




g e (x)e ltx dx = $ (e + it) - $„-(-e + it) 



The strategy of Rudin is to use Cauchy's Theorem to replace $o and <& n by $ 7r /2 
if t < —A and by <&- n /2 for t > A. The problem here is that the function g is not 
entire. We need to avoid the branch cuts to apply Cauchy's Theorem so we define, 
for rj > : 

rg = [0;k-rj] U [k + m +oo)U{k + V e te /e e [0;tt]} and = {s/ - s e T n } 
Then g is analytic on Tq and on T^. 

Let us also denote by &a{w) = J T v g(z)e~ wz dz for a = and ct = tt. The interval 
[A; — rj; k + rj] of T has been replaced by a half-circle to bypass the branch point k 
of the square root so that Rudin's arguments to use Cauchy's Theorem now hold : 
$2(e + it) and $^(e + it) both coincide with ^ 7r / 2 (e + it) if t < -A. (Note that the 
domain of $ 7r / 2 is n„./ 2 = {^/Sz < —A}.) 



Lemma 1 

It holds lim^o $o(w) = Jr g(z)e~ wz dz and lim^ = Jrv g(z)e' wz dz. 

Proof. 

The function c/ is continuous on both half-circles 7^(/c) = {/c + r\e %d jQ G [0; tt]} and 
7»,(-fc) = {~k + ^/fl G [0; tt]} so 



/ g{z)e~ wz dz 



< max \g(z)e~ wz \ I • L(^ v (k)) 



where L(7^(/c)) is the length of the path 7 r? (/c) i.e. n-rj. Now this path of integration 
is compact and the integrand is continuous so its maximum is attained for some zq. 
The half-circle tends to the point k as 77 tends to so zq tends to k and the maximum 
value tends to the value of the function at k which is zero. Then / 7 ^ g(z)e~ wz dz 
tends to zero as rj tends to zero. The same result holds for the path ^y v (—k). And 
the limits of the lemma are a consequence. ■ 
Thus, for any t < —A : 

J IR g e (x)e- itx dx = $ (e + it) - $^(-e + it) 

= lim T; ^o $o(e + it) - lim^o $Z(~ e + **) 
= $ ff/2 (e + it) - $ w/2 (-e + #) 

It follows lim e ^ IjJl9e(x)e~ ltx dx = 0. ■ 

3.3 The main theorem : causality for the nonlinear problem 

The aim of this last section is to prove that the propagation is causal i.e. that a 
compactly supported initial signal with vanishing initial velocity is still compactly 
supported as time goes by. It also means that the propagation of the wave front is 
not faster than light. Deutch and Low assert in [16] that for the linearized problem 
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with a potential barrier, the solution is causally related to the initial condition 
for propagation to the right. Their idea is that the violation of causality is only 
apparent since the energy contained in the forward tail of the initial function (a 
Gaussian centered one) is redistributed. In their situation, the initial condition on 
u t does not vanish but is such that their initial wave packet is moving to the right. 
Our conditions are not so restrictive and our solution is the superposition of right 
and left-moving waves. Our initial condition is supposed to be compactly supported 
which is necessary to prove causality properly. The following theorem states that 
causality still holds when adding a nonlinearity. Our proof is an adaptation of Reed 
and Simon's Theorem X. 76a in [28]. 

Theorem 7 

Assume that f is such that (/i,/2) belongs to D{A) and that the support of f\ is a 
compact set included in (0; R) and f 2 = 0. Then the solution of Problem (NLP), 
given by any of the theorems of the latter section, has the property that the support 
of u(t; •) is a subset of 

C(S, t) = {x G IR/dist{x, S) < ct} 

Proof. 



Lemma 2 

The reflection and transmission coefficients Cr(uj) = (— ° 2 Ql ) and 

\ujc + v c ^ — a 2 + a i J 

T 1 ' 2 {uj) = ( -, 2uJC = ) both define bounded functions of uo in the upper half- 

plane {uj/Im{uj) > 0}. 

Proof. 

Multiplying both coefficients by ujc — y/c 2 u 2 — 02 + a± leads to 

(luc — \Jc 2 ijj 2 — 02 + a\) , „ 2uc (ujc — Jc 2 uj 2 — a2 + eta) 

C R (u) = ± >- and T 1 ' 2 H = ^ >- 

a 2 — ai a 2 — a\ 

The first expression to study is : ujc — \J(?uj 2 — a 2 + a[. Since it is holomorphic in 
the upper half-plane, it is bounded provided that its modulus tends to a finite limit 
when \u\ tends to +00. Denoting by r and 9 the modulus and argument of u, it 
holds : \uc — v 1 c 2 uj 2 — a 2 + a\\ 2 = 

(cr cos(^) - M(r) cos(A(6))) 2 + (cr sin(^) - M(r) sin(A(^))) 2 
with M(r, 9) := ((c 2 r 2 cos 2 (2^) - a 2 + a x ) 2 + c 4 r 4 sin 2 (2^)) 1/4 
and A{r, 6):=\ arctan f /^y) \ 

\ ' / 2 \c 2 r 2 cos(29)—a2+ai J 



Now an easy computation gives M(r, 6) := cr 



(ai - a 2 ) cos(26») / 1 
1 H ho — 

2c 2 r 2 \r 2 

(a 2 — ai) cos(26') /l 



and 



A(r, 9) = 9 so that luc — \/c 2 uj 2 — a 2 + aA = h o 

2cr \r 

Thus lim^^ C R (re ie ) = and lim^^ T li2 (re ifl ) = cos(2#). ■ 

Note that if u is real, the reflection coefficient tends to and the transmission one 
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to 1 as the frequency u tends to infinity which is a reasonable and unsurprising 
behaviour. 



Let us first consider U\ : 

The proof of Reed and Simon (cf. [28], p. 309) is adaptable. Only the first part of 
the proof concerning the linear part has to be changed in the following way. The 
solution reads in terms of Fourier transforms : 



u±(t,x) = F^x [cos(Voj + cWt)?^ 



COS ( di + C 2 U 2 t) 



' UJC — \J C 2 UJ 2 — 0,2 + 0,1 



y uc + V c 2 uj 2 — a 2 + ai 
Let us check the assumptions of Theorem 6 for g defined by 



</(*) = 



cos (yj ai + c 2 z 2 t) 



cz — y/c 2 z 2 — a 2 + di 
cz + \/c 2 z 2 — a 2 + a{ 



Ffi{z) 



First of all, since /i is an L 2 -function with compact support contained in (0;R), the 
Fourier transform of f\ is an entire analytic function and there exists a real constant 
C\ such that 

\TlM\<CieW R 
(cf. [31], introduction of Chapter 19 for this classical result). 

Further since the roots drop out in the power series, cos(V a i + c 2 u 2 t) is also entire 
and satisfy, for some real constant C 2 



cos 



< C 2 e 



\u\ct 



Now, since the modulus of the reflection coefficient is bounded (cf. Lemma 2), there 
is a real constant C such that 

\g{uj)\ < Ce^ R+ ^ 
i.e. the estimate of Theorem 6 is satisfied with A := R + ct. 

Furthermore, since f\ is an L 2 -function, then its Fourier transform lies in L 2 (IR) as 
well. And since the reflection coefficient is bounded on IR as well as the cosinus, g 
belongs to L 2 (IR). 

Since the reflection coefficient is an analytic function in the upper half-plane only 
and it has two branch cuts : j(±k) with k = \fd 2 — a±, the function G defined on IR 
by G(x) = jjfig(uo)e~ lulx doj is such that the support of G is a subset of (—A; +oo) 
by Theorem 6. Then that of x i— > G(— x) is included in (— oo; A) i.e. (— oo; R + ct). 
Since Ui(t,-) is the restriction to IR + of this function, it is compactly supported 
in (0;R + ct) i.e. its support is in C((—R; R), ct). The rest of the proof remains 
unchanged. 



As for u 2 , it is rewritten as : u 2 (t,x) = 
■5^ f_ + ~ cos(Vai + cWty^ei-"-^ 



c 2 u) 2 —a2+ai 



)-( 



2wc 



Tji (uj)du 



i ^ — UJ— i \J C 2 U) 2 — 0,2+0,1^ X 



Now, since the transmission coefficient and the additional factor e 
are bounded (cf. Lemma 2), the same arguments as before give the compactness of 
the support of u 2 as regards with the space variable x. m 
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